Bulk viscosity and the conformal anomaly in the pion gas 
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We calculate the bulk viscosity of the massive pion gas within Unitarized Chiral Perturbation 
Theory. We obtain a low temperature peak arising from explicit conformal breaking due to the pion 
mass and another peak near the critical temperature, dominated by the conformal anomaly through 
gluon condensate terms. The correlation between bulk viscosity and conformal breaking supports a 
recent QCD proposal. We discuss the role of resonances, heavier states and large- A^c counting. 
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The matter produced after thermahzation in relativis- 
tic heavy ion coUisions behaves nearly as a perfect fluid 
[T]. Deviations are seen mainly in elliptic flow and can 
be reasonably explained with a small shear viscosity over 
entropy density ratio r]/s < 0.5 [5], whereas bulk viscos- 
ity ( is generically assumed to be negligible. However, 
it has been recently proposed [3] that C, might be large 
near the QCD phase transition. If C/s is comparable to 
ry/s near the critical point (where indeed the latter is ex- 
pected to have a minimum) interesting physical possibili- 
ties arise such as radial flow suppression, modifications of 
the hadronization mechanism [3] or clustering at freeze- 
out [4j. The argument of [3] is that, following the QCD 
sum rules in [S], one can relate with the trace anomaly: 
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with {9)t = {TI^)t = e — 3P, T/,' the energy- momentum 
tensor, e the energy density, P the pressure and eq — 
{9)q/A in vacuum. To derive ([ij, a particular ansatz 
has been used for p(w), the {09) spectral function at 
zero spatial momentum, with (p/ti;)(0) = 9(^/7r and 
QC,LOo = 2 {p/uj)duj. Eq. (nj) implies then a large 
bulk viscosity near the QCD transition, from the {9)t 
peak observed in the lattice [6^ , more or less pronounced 
depending on the transition order [3^. However, this ar- 
gument has been recently criticized on the basis of the 
(p/lu) convergence and parametric dependence with 
the QCD coupling constant [7J. On the other hand, esti- 
mates of C from lattice data show that ujS{uj) terms and 
large-o; non-thermal contributions have to be properly 
accounted for in spectral functions [8, 9J. 

It is therefore of great importance to study QCD 
regimes where one can rely on analytic calculations, in 
order to clarify the validity of the above proposal without 
appealing directly to lattice data. In the weak couphng 
regime, valid for very high temperatures, (^/r] has been 
found to be parametrically small ^lOj . Another regime 
where one can perform analytic calculations is low-energy 
QCD, where the system consists primarily of a meson gas 
and, for low temperatures, one can rely on Chiral Per- 
turbation Theory (ChPT) [11]. In this regime, we have 
recently shown |12l I13| , within Linear Response Theory 



(Kubo's formulae), that the usual ChPT power count- 
ing must be extended to account for 1 / Fp contributions 
arising in transport coefficients. Here, Fp is the thermal 
width of a pion with three-momentum p, in which the tttt 
total elastic cross section enters linearly in the dilute gas 
regime |14j . Performing the power counting, which in- 
cludes a detailed analysis of ladder-type diagrams consid- 
ered in [15], the leading-order ChPT contribution comes 
from a one-loop meson diagram with Fp internal 
lines. An essential point is to include unitarity correc- 
tions in Fp to describe correctly the temperature behavior 
as the system approaches chiral restoration. We neglect 
inelastic 2tt <-!■ An reactions restoring particle number 
equilibrium, which are suppressed in our counting and 
yield chemical relaxation times about ten times larger 
than the plasma lifetime [TB]- Thus, our bulk viscos- 
ity is meaningful for the pion gas formed in heavy ion 
collisions, which conserves approximately pion number 
between chemical and thermal freeze-out, as confirmed 
by particle spectra data analyses with a pion chemical 
potential |17| . If C is defined in complete chemical equi- 
librium, then particle-changing processes dominate [15j . 
The dominance of elastic processes for ^ in the pion gas 
holds also in kinetic theory [TH [T9j EO]. With our ap- 
proach we have also obtained rj/s developing a minimum 
compatible with AdS/CFT bounds with values in good 
agreement with kinetic theory |19l I21| and phenomeno- 
logical estimates on elliptic flow. This is the theoretical 
basis of the present work, where we will analyze within 
ChPT the correlation between bulk viscosity and the con- 
formal anomaly in the pion gas regime, studying the ori- 
gin of the different contributions to conformal breaking 
for physical massive pions. Thus, we start with Kubo's 
formula: 
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where the modified pressure operator V = — T?/3 — c^Tqo, 
the squared speed of sound — dP/de — s/cy, s — 
dP/dT and the specific heat c„ = de/dT = Tds/dT. 
We follow the conventions of [22], where C is defined as 
the change in the pressure produced by a gradient in the 
fiow velocity, relative to equilibrium. This leads to the 
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correlator in ([2]), which is the adequate one to be used 
within perturbation theory [Tj [10]. In lattice analyses, 
one works with the Lorentz invariant 6 instead. In our 
approach these two correlators are not equivalent, since 
the leading order in l/Fp for perturbative Too commuta- 
tors does not vanish for zero spatial momentum. As we 
shall see, sticking to the original definition ([2| leads nat- 
urally to the expected conformal properties and asymp- 
totic behavior of the bulk viscosity. Following [T^], we 
calculate then the spectral function (V commutator) in 
([2]) in the imaginary-time formalism, picking up the dom- 
inant contribution in 1 / Fp (pinching pole) of the analyt- 
ically continued retarded correlator. That term is purely 
imaginary and gives the dominant effect in the spectral 
function at zero momentum and small energy. Thus, to 
leading order: 

^ ^ 

(3) 

with n-Q{x) = l/(exp(a;/T) — 1) the Bose-Einstein distri- 
bution function, Ep = yjp^ + and where the leading 
O(p^) order in T^'^ has been retained in the vertex. Now, 
we get cl in ^ from P calculated up to 0{T^) in [TT] . 
In Fig|l]we see that to O(T^), both the specific heat and 
the speed of sound increase monotonically, approach- 
ing the ultra-relativistic limit of 1/3 corresponding to a 
gas of free massless pions. To that order, since the dis- 
tribution function is peaked around p ^ T for T 3> tjIt^, 
we see that ([s]) vanishes asymptotically for large tem- 
peratures, as expected for conformally invariant systems 
[M ESI mi m HI- In fact, from ^ we get for mass- 
less pions (chiral limit) C — 15(1/3 — 0^)^77, consistently 
with [2,3 and paramctrically with high-T QCD [10 . The 
crucial point here is that taking one more order in the 
pressure Cy grows, reaching a maximum at about ~ 
220 MeV. The speed of sound attains then a minimum 
at Tc which will alter the behavior of C{T). This is the 
critical behavior of a 0(4)-like crossover, as expected for 
two massive flavors at zero chemical potential. A physi- 
cal interpretation is that, although temperature tends to 
erase mass scales, chiral interactions are enhanced and 
produce in the critical region a significant, nonpertur- 
bative, conformal breaking reflected in ^ 1/3. Note 
that, although in the massive case Tc is near the chiral 
restoration temperature where the order parameter 
{qq)T vanishes [TT], in the chiral limit ~ 170 MeV, 
while Tc is almost unchanged. 

We plot our result for the bulk viscosity in Figj2] The 
effect of including the O(T^) in effectively produces a 
peak around Tc, not present to 0{T^) [13]. The speed 
of sound is not the only relevant effect yielding a siz- 
able peak: unitarization of the cross section entering Fp 
[T^ is also crucial to O(T^). Considering unitarized par- 
tial waves for tttt scattering (ChPT is only perturbatively 
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FIG. 1: Specific heat (left) and speed of sound squared (right) 
for the pion gas. The red dashed line is the 0{T^) calculation, 
and the continuous blue line the 0{T^) one. The green dotted 
line is the ultra-relativistic limit — 1/3. The dashed-dotted 
blue line is the chiral limit result to C'(r**). 

unitary) improves the high energy behavior (and there- 
fore the high temperature one) and generates dynami- 
cally the /o(600) and p(770) resonance poles. Consis- 
tently, we have chosen the values of the low-energy con- 
stants k entering pion scattering (they can be found in 
[2^) so that the mass and width of the p are at their 
physical values for T = 0. As we discuss below, the k de- 
pendence is crucial in the present analysis. In the chiral 
limit, the transition peak is almost unchanged and so is 
Tc, unlike T^, which indicates that chiral restoration is 
not the main source of this effect. Our massless results 
are in reasonable agreement with a recent kinetic the- 
ory analysis |20j. We also obtain a low-T peak, which 
disappears in the chiral limit. In our regime and for 

momenta p = 0{\/rr^/T) and taking the leading order 
for Fp [12J and ~ T/m^ -|- . . . [U], eq.(|3]) becomes: 

C(T)~13.3^^^ for r«m^, (4) 

where is the pion decay constant. The above behavior 
is consistent with nonrelativistic kinetic theory |18j where 
C and 77 are expected to be comparable at low T. Thus, 
C(T) increases for very low T and has to decrease at some 
point to match the asymptotic vanishing behavior, thus 
explaining the low-T maximum. 

Let us now evaluate conformal-breaking contributions 
for the pion gas. First, it is instructive to recall the QCD 
result for the trace anomaly [53]: 

i^i^)QCD = ^f^G^l^Cr + (1 + lm{g))qMq, (5) 

where the renormalization group functions are, pertur- 
batively, I3{g) = 0(5^), 7m(5) = 0{g^). The flrst term is 
the conformal anomaly proportional to the gluon conden- 
sate. The second one comes from the explicit breaking 
in the QCD lagrangian, M being the quark mass matrix. 
For the pion gas, using the thermodynamic identity 
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FIG. 2: Bulk viscosity of the pion gas. The full blue line is the 
unitarized result with to 0{T^) and the dashed-dotted blue 
one is the same calculation in the chiral limit. The dashed 
green line is the nonunitarized result at the same order. The 
dotted red line is unitarized with to 0{T^) and lies very 
close to the nonunitarized curve, which is not displayed. 

we represent in Fig|3]the trace anomaly to different orders 
in the pressure, as well as the T-function appearing in 
the r.h.s. of Q. We observe clearly the same two-peak 
structure as the bulk viscosity, with similar features. 

The low-T peak disappears in the chiral limit. Its 
contribution comes then from explicit conformal break- 
ing. Calculating only the first nonvanishing order in 
ChPT, either using ^ or evaluating directly the energy- 
momentum correlators, we get: 

{0)t-{0)o - 3mlg,{m^,T) + Oif-^) 

^ 2m,i{qq)T-{qq)o) + Oif-^), (7) 

where niq — rriu — rud and we formally account for differ- 
ent chiral orders by their power. The function gi is the 
thermal correction to the free pion propagator G{x — 0) 
Comparing with the QCD expression ([5| the factor 
of two in ([t]) for the quark condensate is perfectly con- 
sistent with the result [26] showing that the quark and 
gluon contributions to the trace anomaly are identical at 
low temperatures. Now, gi{T)/T'^ has a maximum at 
T ~ 2to^/5 ~ 60 MeV, which is the low-T peak in Fig(3] 
and the source for the first peak of the bulk viscosity. 

The transition peak only shows up at 0{T^) and sur- 
vives in the chiral limit, where its origin is purely anoma- 
lous. It comes from ChPT interactions involving dimen- 
sionful couplings, like and is therefore suppressed at 
low temperatures |27j . For massive pions, the value of 
the peak and its position are almost unchanged with re- 
spect to the chiral limit, the difference being even smaller 
than the simple extrapolation of the quark condensate 
contribution in Q with {qq)T to 0(T*), which repre- 
sents around a 10% correction in the critical region. The 
fermion contribution is also subdominant in lattice anal- 
yses [6]. These results show again that the nature of 
this effect is not likely to be related to chiral symmetry 
restoration but rather to other QCD critical effects like 
deconfinement. The correlation with the bulk viscosity 
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FIG. 3: Thermal expectation value of the trace anomaly for 
the pion gas. The dashed-dotted green and continuous blue 
lines are, respectively, the 0{T^) and 0{T^) results. The 
dashed red line corresponds to the 0{T^) result for massless 
pions (the 0{T^) order vanishes for — 0). The dotted 
blue line is T^^^^i^. 



is again clear. In fact, in the chiral limit the function be- 
tween brackets in ([l]) and 1h{c1 — 1/3)^ = Q/rj have their 
maximum at the same = e~^/^ Ap with Ap given in [TP 
in terms of li + 4^2. We recall that in order to establish 
the possible correlations between the conformal anomaly 
and the bulk viscosity, we have used the same set of U 
in both figures. For those unitarized values, — 220 
MeV. Using perturbative values, for instance those given 
in [TT] fixed to reproduce pion scattering lengths, the 
critical peak is about three times smaller and Tc — 148 
MeV, while varies only about 10 MeV from one set 
to another. We get exactly the same drastic reduction of 
the critical peak and shift of Tc in the bulk viscosity. The 
presence of resonances is then crucial to yield a sizable ef- 
fect in the transition peak, whose dominant contribution 
comes from the gluon condensate. 

Regarding the u}o{T) function defined through ([TJ, 
in the chiral limit it grows linearly with T, reaching 
tJo ~ 400 MeV at the transition. In the massive case, 
taking |eo| — /^'tiJ, the ChPT lowest order, we get 
ujo{Tc) ^ 1 GeV, almost constant from T - 150 MeV 
onwards. These values are in reasonable agreement with 
the estimates in 3 . On the other hand, from Q we get 
LOo{T) ~ OASmV^if^VT) for T -> 0+. 

The numerical values of the trace anomaly in Figj3] 
are not far from the lattice values [5] for low T, but 
they are about a factor of 10 smaller near Tc. The in- 
creasing of degrees of freedom due to heavier states, not 
included in our approach, is clearly important in that 
region. For instance, the 0(T®) pressure in the chiral 
limit is proportional to N'j{N'j — 1) [571 so that chang- 
ing from two to three flavors, which are not Boltzmann 
suppressed near the transition, increases significatively 
the anomaly. In addition, using a simple free Hadron 
Resonance Gas approach [3], the 7r,p,cr contributions to 
the anomaly amount only to a 5% of all baryon and 
meson states up to 2.5 GeV. In fact, although we get 
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C/s ~ 0.02 at the transition peak, still much smaller than 
ry/s ~ 0.25, we would get a larger value if we assume that 
the introduction of heavier states increase the anomaly, 
and that implies an increase of the transition strength 
and a strong reduction of [3]. As an indication, set- 
ting = in ([3| we get C/s ~ 1 at T^. 

We have seen that it is crucial to include correctly 
the effect of the p resonance. On the other hand, the 
/o(600)/cr is expected to be related to chiral restoration. 
Regarding bulk viscosity, it has been suggested in [2H], 
within mean field theory, that any dynamic scalar field a 
should contribute to ^ cx Fcr/m?^, which may be large near 
the critical region by mass reduction, for instance in the 
Linear Sigma Model (LSM) context. Within unitarized 
ChPT, the dynamically generated /o(600) pole under- 
goes a significant mass reduction towards 2m7r governed 
by chiral restoration, remaining a broad state with sizable 
width near the transition [Mj. Interestingly, from [23], we 
find that Fa/ml. has a peak at T ~ 180 MeV, where the 
pole mass reaches threshold. For higher T the width still 
decreases (by phase space reduction) while the mass re- 
mains close to threshold. This critical value is very close 
to the one obtained in ^S] for the LSM assuming a T- 
independent width. However, as discussed above, these 
chiral restoration effects are likely to be subdominant. 

The large-iVc limit is also revealing. The counting of 
the li can be extracted from the Li {Nj — 3) [53] while 
/2 = 0{Nc). We get Fp ~ 0{N-^) and, in the chiral 
limit, {9)t - 0{N-^) - (c2 - 1/3) so that C - C(l) 
and (/rj ~ 0{N^^), parametrically suppressed as ex- 
pected. Now, taking into account the critical behavior, 
Te ~ 0(e^=) and (6I)t, ^ O^e'^^ /N^). This large depen- 
dence is another indication of the dominance of confine- 
ment over chiral restoration, comparing with the chiral 
= 0{Nc). Also, {6) cx L3, which in large-TVc in- 
cludes a term proportional to the gluon condensate [30] . 
Comparing with the QCD expressions in [lOj . we agree 
except for the overall 0{N^) constants in the pressure 
which count the degrees of freedom. For massive pions, 
the above chiral limit scaling is only reached asymptoti- 
cally for large F, while for any T we get C/rj ~ 0{1) ^ 
{0)t - {0)o with C 0{N^), compatible with Q. 

Summarizing, we have shown, within unitarized ChPT, 
that the massive pion gas develops a strong correlation 
between bulk viscosity and the conformal anomaly. Both 
quantities show a low-temperature peak coming from 
mass conformal breaking and another one at the criti- 
cal temperature remaining in the chiral limit and mainly 
dominated by gluon condensate contributions not re- 
lated to chiral restoration. The dynamically generated 
light resonances are essential to obtain sizable effects at 
the transition. Different estimates indicate that heavier 
states could yield a larger bulk viscosity near the transi- 
tion, leading to observable effects in heavy ion collisions. 
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